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IGEC2: A 17-month search for gravitational wave bursts in 2005-2007 
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We present here the results of a 515 day search for short burst of gravitational waves by the IGEC2 
observatory. This network included 4 cryogenic resonant-bar detectors: AURIGA, EXPLORER 
and NAUTILUS in Europe, and ALLEGRO in America. These results cover the time period from 
Nov 6 2005 until Apr 15 2007, partly overlapping the first long term observations by the LIGO 
interferometeric detectors. The observatory operated with high duty cycle, namely 57% for 4-fold 
coincident observations, and 94% for 3-fold observations. The sensitivity was the best ever obtained 
by a bar network: we could detect impulsive events with a burst strain amplitude h rss ~ 1 ■ 10~ 19 
with an efficiency ^,50%. The network data analysis was based on time coincidence searches over at 
least three detectors, used a blind search technique and was tuned to achieve a false alarm rate of 
1/century. When the blinding was removed, no gravitational wave candidate was found. 

PACS numbers: 04.80.Nn, 95.30.Sf, 95.85. Sz 



I. INTRODUCTION 

In the last decade, while interferometric detectors were 
approaching, through a series of commissioning runs, 
their target sensitivity, resonant bar detectors contin- 
ued to reliably observe the cosmos for very long periods 
of time, looking for unmodeled impulsive gravitational 
waveforms (GW) from galactic sources. 
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For the purpose of coordinating this search effort, 
the four cryogenic, resonant mass detectors (plus, in 
an early stage, the Australian antenna Niobe) joined 
forces in a collaborative agreement called IGEC (Inter- 
national Gravitational Event Collaboration) that mainly 
consisted of a protocol for data exchange and analysis. 
The IGEC collaboration analyzed almost 4 years of data 
from 5 antennas This protocol was later revised 

and renamed IGEC2: among other changes, it included 
coordination on scheduling routine maintenance opera- 
tions (refills of cryogenic fluids) in order to ensure max- 
imal time coverage with at least 3 of the 4 detectors. A 



2 



first period of 6 months of data was analyzed and pub- 
lished under the new protocol [3J: it yielded a negative 
result (no GW candidates) but it was instrumental in 
setting up and testing network analysis procedures that 
were the starting point from which the present analysis 
has evolved. 

The detectors taking part in the IGEC2 network 
are: ALLEGRO, located at Louisiana State University 
(Louisiana-USA) and operated by the local ALLEGRO 
group, AURIGA, located in the Legnaro National Lab- 
oratories (Padova-Italy) of INFN and operated by the 
AURIGA collaboration, EXPLORER, located at CERN 
(Geneva-Switzerland) and NAUTILUS, located in the 
Frascati National Laboratories (Frascati-Italy) of INFN. 
Both latter detectors are operated by the ROG collabora- 
tion. AURIGA, EXPLORER and NAUTILUS all benefit 
from support from INFN, while ALLEGRO was funded 
by NSF. 

This network of four detectors, all operating with high 
duty cycle and large overlap time, is effective in reject- 
ing spurious candidate events caused by transient local 
disturbances or by intrinsic detector noise. The sensitiv- 
ity of these antennas has been superseded by the much 
superior performance of large interferometers like LIGO 
and VIRGO. Therefore, IGEC2 upper limits on the GW 
flux on Earth are no longer astrophysically significant so 
in this paper we target the possible detection of a rare, 
impulsive event with long term observations. 

In the previous runs of this network, either the sensi- 
tivity of the antennas was lower pQ [2] , or the observation 
lasted for a much shorter stretch of time [3J- The LIGO 
observatory has also carried out searches for bursts, both 
with its three interferometers i4j|5jj6] and in a coinci- 
dence run with GEO [7] . 

The first published results |9] of the LIGO S5 run 
report a relevant advance in these searches, yielding a sig- 
nificant lowering in the previously published upper limits. 
That work refers to the period November 2005 through 
November 2006, covered with a ~ 70% duty cycle. The 
same period is covered by the present work, clearly with 
a lower sensitivity, but with a much larger duty cycle 
(~ 94%). For this reason it can be used in conjunction 
with this early phase of S5, when the most common pe- 
riods were stretches of time with only one interferometer 
taking data. 

Several astrophysical processes can generate gravita- 
tional radiation in the sensitive bandwidth of the IGEC2 
detectors (all centered around 900 Hz): among these, 
stellar core collapse [TU], final phase of inspiralling and 
ring down in merger of compact binary systems, as well 
as rotational instabilities and quasi-normal mode oscil- 
lations of relativistic stars [11] [12]. The amplitude of 
these events, according to current estimates [T3] are such 
that an event in our galaxy (say within a few kpc) can be 
detected by the resonant antennas. In view of the large 
uncertainty in the modeling of these waveforms and of 
the restricted bandwidth of our detectors, we search our 
data for unmodclcd, featureless short bursts: as in the 



previous search [3J, we aimed at signals with typical du- 
rations up to a few tens of milliseconds (see the analysis 
of sect III I, e.g. damped sinusoids, ring down, or gaussian 
bursts. 

The main advances reported here refer to a detailed 
analysis on the choices of the coincidence window and of 
the thresholds strategy, carried out through an extensive 
use of software injections of pulses in the data. In addi- 
tion, we analyzed a substantial stretch of time of fourfold 
operation. This paper describes the results of 17 months 
of observation, from Nov. 16 th 2005 to Apr. U th 2007, 
containing the longest reported period of fourfold coinci- 
dence observation, 293.5 days out of 515. Indeed, in the 
first IGEC search, published in 2003 0, only 26 days of 
fourfold observation were collected at a much lower sensi- 
tivity, in a four year span and using five antennas. In the 
2007 paper[3J we only analyzed threefold coincidences. 

The goal of this observation campaign, just as in that 
reported in [3] was to look for coincidence events with 
a background of accidentals of 1 per century. This low 
level of false alarm was achieved with a suitable choice of 
the thresholds applied to the candidate event lists pro- 
duced by each group, as discussed in detail in sect |III| Wc 
chose to undertake five searches in parallel: one four-fold 
coincidence search and four triple coincidence searches, 
a priori assigning to each a false alarm rate (FAR) of 
0.2 / 'century. Once these five searches are combined in a 
logical OR, the desired level of FAR is achieved. Two-fold 
coincidences were not analyzed, since in order to achieve 
a reasonable value of the FAR, we should have had to 
raise the thresholds so much that the efficiency would 
have been even lower than for the threefolds. 

In the following section we review some features of the 
network, of its detectors and of the procedure used to 
produce the exchanged data. 

In sect |III| we describe the study carried out to eval- 
uate the background and set up the corresponding cuts 
in the data that guarantee the desired FAR. Finally the 
results of the search are presented and discussed in the 
last section. 



II. THE IGEC2 OBSERVATORY: THE FOUR 
DETECTORS AND THEIR DATA 

The four cryogenic resonant GW detectors that are 
part of the IGEC2 network are reliable and stable ma- 
chines. Their long term noise performance is rather sta- 
tionary. Therefore, both the general observatory descrip- 
tion and the noise characterization of the antennas are 
substantially unchanged with respect to that described in 
[3J. While detailed description of the equipment can be 
found in [HITS] for AURIGA, in [H[T7] for EXPLORER 
and NAUTILUS and in [TH] for the ALLEGRO detector, 
respectively, here we will just recall the common features 
of the antennas. 

These four detectors are very similar in their experi- 
mental set up. They all consist of a 3 m long Aluminum 
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bar, suspended in vacuum in a cryogenic environment. 
The quantity we monitor and analyze, in search of an 
event generated by GW, is the amplitude of vibration 
of the first longitudinal elastic mode of the cylinder, res- 
onating around 900 Hz. This signal is intrinsically narrow 
band, the bandwidth being set by the interplay between 
thermal (resonant) noise and amplifier (wide-band) noise 
and reaching at best 10% of the operating frequency. For 
this reason, the signal we can extract from the antenna 
vibration can be related to the Fourier amplitude H(u>) 
at the resonant frequency of an hypothetical short GW 
signal, rather than to the wave amplitude h(t) itself. 

The spectral sensitivity curves of the four antennas are 
shown in figjl] Different spectral shapes in the various 
detectors arise from particular choices of the antenna- 
readout coupling. Indeed, the resonant transducer that 
is used to convert the bar vibration into an electric sig- 
nal (a light mass mechanical oscillator, whose vibrations 
modulate an electric or magnetic field), is different in 
the four detectors. While ALLEGRO has a supercon- 
ducting, persistent-current, inductance-modulation de- 
vice, AURIGA, NAUTILUS and EXPLORER rely on 
a capacitive, constant charge biased transducer. In all 
detectors the first stage amplifier is a d.c. SQUID, an 
ultra-low noise device. AURIGA indeed uses a double 
d.c. SQUID with the matching LC circuit tuned to the 
mechanical antenna frequency (resulting in a three mode 
resonant system). On the other hand, EXPLORER and 
NAUTILUS, while using a loosely tuned [T7] matching 
circuit, employ a very small gap capacitor (8 — 12/im) 
to achieve a high electromechanical coupling. In ALLE- 
GRO, due to the low source impedance of the inductive 
transducer, no matching transformer is needed. Both 
NAUTILUS and AURIGA are equipped with dilution re- 
frigerators that allow them to operate at ultralow tem- 
perature. However, because of reliability and duty cycle 
considerations, they were in this case operated at about 
3 - 4 K, like EXPLORER and ALLEGRO. With respect 
to the spectral sensitivity curves shown in [3], the Sh(f) 
of Explorer now has a more symmetric shape. This is due 
to a small change in the bias voltage of the resonant ca- 
pacitive transducer that was implemented in April 2006. 

All four detectors require periodic down-time for cryo- 
genic maintenance, typically one or two days per month 
for refilling cryogenic fluids. A coordination effort was 
made to undertake maintenance duties on different days 
for each apparatus, in order to maximize the time with 
at least three detectors simultaneously in operation. As 
a result, we achieved the observation times described in 
table HJ that shows that we had three or more antennas 
operating for 94% of the entire 515 day period with data 
quality suitable for science goals. 

This long period of observation and the high duty cycle 
of all detectors guarantee, despite the periodic mainte- 
nance stops, an adequate coverage of all hours in the 
diurnal time, as shown in fig(2] where the percentage 
of time of the fourfold mode of operation is plotted vs 
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FIG. 1: Typical strain noise spectral densities (single-sided) 
of IGEC2 detectors. All detectors are sensitive in a region 
around 900 Hz, and the minimum level of noise is compara- 
ble in all spectra. The wider bandwidth of AURIGA includes 
the bandwidths of the other detectors. The noise curve of 
Explorer (blue on-line) refers to operation after the 2006 re- 
tuning. 



Configuration 


Time of operation 
(days) 


detector 





1 detectors 


1.6 


2 detectors 


31.0 


3 detectors 


188.8 


4 detectors 


293.5 



TABLE I: Multiplicity of observing antennas during the ana- 
lyzed period. 

the hour of the day. While the sidereal time distribu- 
tion has small fluctuations around the expected value of 
4.17%, the solar time distribution exhibits larger oscil- 
lations with a dip in coincidence with the start of local 
morning activities in Italy and Geneva. 

As in previous joint searches [3J, each group is respon- 
sible for the production and the calibration of the data 
of its own detector. This includes independent design 
and operation of the antenna hardware, data sampling 
and filtering and extraction of the candidate events, by 
proper choice of the amplitude thresholds. All groups fil- 
ter the data by applying adaptive linear filters matched 
to (5-like signals. We carried out calibrations of these 
filters based both on mechanical excitations with short 
pulses and on software injections of various impulsive 
waveforms. These calibrations show that the (5-filtered 
data preserve a good efficiency also to longer signals, such 
as damped sinusoids with decay times r < 30 ms, as it 
will be discussed in more details in sec.IIIA. 

A candidate event is identified by detecting a lo- 
cal maximum in the absolute value of the filtered data 
stream. The occurrence time of the maximum is the es- 
timate of the arrival time of the burst and its amplitude 
is the estimate of the Fourier amplitude H of the GW 
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LEGRO ceased operations. 



Detector 


Noise 

(Hz- 1 ) 


Threshold 


Number 
of events 


ALLEGRO 


2.7 • 10~ 22 


zr = i.i- w-' 21 h z - l 


1,472,517 


AURIGA 


1.5- 10~ 22 


SNR=4.5 


585,968 


EXPLORER 


4.1 • 10~ 22 


SNR=4 


1,193,830 


NAUTILUS 


3.5- KT 22 


SNR=4 


1,400,882 



TABLE II: The average noise level, the value of threshold 
chosen by the responsible group, and the resulting number of 
exchanged events, for the four detectors. 
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FIG. 2: Distribution of the fourfold observation time in the 
hours of the day: sidereal (solid line, red) and solar (dotted 
line, blue). 

waveform h(t), assumed to be a delta. In order to re- 
move short time clustering in the events of each detector, 
we introduced a dead time in the process of selection of 
the events. The dead time represents the minimum time 
allowed between two different events, and its value was 
chosen 1 s for EXPLORER, NAUTILUS and AURIGA 
and 0.2 s for ALLEGRO. 

Thus, AURIGA, EXPLORER and NAUTILUS adap- 
tively set their thresholds, depending on the local (in 
time) noise level. However, due to a different data pro- 
duction process, ALLEGRO, whose data where not used 
in the 2005 search, produced its own events with a fixed 
H threshold. 

The SNR (or H for ALLEGRO) threshold value and 
the resulting number of event candidates, as well as the 
average noise level of each detector are listed in table 
pi} Note that, although independently chosen by each 
group, the applied thresholds yield a comparable num- 
ber of events (1 to 3 thousand per day) from all detec- 
tors. The ampllitude distribution of all exchanged events 
is shown in fig. [3] 

Each group exchanged the events list after adding a 
time offset, kept confidential to the other groups. Back- 
ground evaluation and the tuning of the time coincidence 
analysis were therefore completed in a blind manner; only 
after agreeing on the thresholds, coincidence windows 
and other choices that contribute the search strategy 
were the confidential shifts revealed (opening the box), 
thus comparing the real candidate event times and find- 
ing the actual coincidences in the data. 

The data exchange and background evaluation was car- 
ried out in two installments: first the event candidates 
for the period Nov 16th 2005 to Dec 31st 2006 (411 days) 
were exchanged and analyzed. The background was eval- 
uated and thresholds were set in order to obtain a FAR 
of 1/ccntury. We call this Data Stretch A. After opening 
the box, 104 more days of data, hence referred to as Data 
Stretch B, became available. As of April 15th, 2007, AL- 



We thought it reasonable to consolidate in this paper 
the entire duration of four-fold operations, but this re- 
quired a separate analysis for stretch B, as that for stretch 
A had already been completed, including the exchange of 
the confidential time shift. In section ITTTl we describe the 
background evaluation carried out separately for the two 
data stretches. 




Amplitude in H {Hz ) 



FIG. 3: Amplitude distribution of the exchanged candidate 
events above the minimal thresholds (AURIGA SNR > 4.5, 
ALLEGRO H > 4 ■ lO -21 /^ -1 , EXPLORER and NAU- 
TILUS SNR > 4.0). 



III. NETWORK DATA ANALYSIS 

We ran five separate searches of coincidences: one four- 
fold, and four triples. Their observation times are well 
overlapped and we combined them in a logical OR, to im- 
prove the overall detection efficiency. In order to achieve 
the desired total FAR of 1 /century, we decided to equally 
partition it into each separate search. 

Once each group has produced its list of candidate 
events, the only parameters that can be adjusted to set 
the desired number of accidentals, i.e. 0.2/century in 
each search, are the coincidence window and a set of sec- 
ond (obviously higher) thresholds. 
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A. The time coincidence window 

In the past [2] [3] , we used approximate analytical mod- 
els to estimate the time response of any detector to delta- 
like signals. 

The time coincidence window was then set, with the 
goal of keeping the efficiency to 95%, by using the esti- 
mated time response within the Tchebichev inequality. 

For this analysis however, we were able to take ad- 
vantage of an extensive methodological study jointly car- 
ried out with the VIRGO group 19J . In this study, we 
applied, to one 2005 day of real data, a large number 
of software injections representing damped sine waves (a 
very convenient and widely used representation of short 
bursts) with various center frequencies (between 866 and 
946 Hz) and various durations (r = 1, 3, 10 and 30 ms): 

h(t) = h sin(27rfot)e- t/T for t>0 (1) 

We summarize here the main results of this analysis, 
relevant for this IGEC2 search: 

• The shortest bursts (r = 1ms) excite the detec- 
tors in a fashion very close to that expected for 
delta pulses, confirming the predictions of analyti- 
cal models, and providing further details and con- 
fidence. 

• The detection efficiency of the antennas remains 
good also for longer burst, provided that they 
carry a sufficient Fourier component in the detec- 
tor bandwidth. If we compare pulses with the same 
h rss , waves with r > 10 ms often produce larger re- 
sponses than the shorter ones, despite the use of the 
delta-matched filter. 

• A large fraction of the timing uncertainty is usually 
due to a systematic offset that has a different de- 
pendence on signal frequency and r in each antenna 
(see figj4|) 

• The width of the timing uncertainty (around the 
offset value) depends on the SNR of the signal and 
only slightly, if at all, on its shape. 

The goal of this study is to set a coincidence window 
as narrow as possible to minimize the background, while 
accommodating all these offsets and fluctuations in or- 
der to keep the detection efficiency high also for non 
delta-like excitations. To this aim, we found it useful 
to investigate the effect of these uncertainties on a pair 
of detectors, rather than on a single one. As an exam- 
ple, we report studies carried out for the EXPLORER- 
NAUTILUS pair. Regarding the systematic time offset, 
fig |4] shows that, as the center frequency or the duration 
of the injected signal is changed, the offset difference can 
vary between -10 and +15 ms. As far as statistical fluc- 
tuations are concerned, fig|5] shows the cumulative dis- 
tribution of the detection time difference in the case of 
short signals. We can see that the first part of the curve 
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FIG. 4: Offsets in the event time differences between Explorer 
and Nautilus for damped sinusoids with r = 1,3,10 and 30 ms. 
This systematic error varies between -15 and +10 ms. 



can be well fitted by a gaussian with a t — 9 ms, while for 
larger time differences there is a longer tail, containing 
about 5 -10% of the events. We can conservatively state 
that 95% of the coincidences is retrieved with a window 
of 25 ms. 

As a consequence of these considerations, we have set 
the coincidence time window to 40 ms, thus including 
the possible offsets deriving from non-5 waveforms, the 
statistical fluctuations and the time of flight (up to ~ 
2 ms for the European antennas). In the coincidences 
involving ALLEGRO data, the window is expanded to 
60 ms, in order to take into account the larger maximum 
(~ 20 ms) time of flight of a signal between ALLEGRO 
and the other detectors. 



B. The choice of thresholds 

In a previous search [3] that focused on triple coin- 
cidences, we explored three different strategies for the 
choice of the thresholds: A) constant, equal SNR in all 
antennas, B) constant SNR in each antenna and higher 
threshold on the more sensitive detector C) constant, 
equal absolute amplitude (H) reconstructed by each an- 
tenna. 

In the present analysis, thanks to the considerations 
detailed below, we were able to focus the analysis on 
one criterion, choosing, since the beginning, strategy B: 
higher SNR threshold for AURIGA, same SNR thresh- 
olds for the other three detectors. This allows, on aver- 
age, a better efficiency to different classes of signal wave- 
forms. 

In fact, in the methodological study cited above [19] . 
we also investigated how different detectors respond to 
the same, non-delta excitation. As an example we show 
in fig|6] the responses of EXPLORER and NAUTILUS 
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FIG. 5: Cumulative number of coincidence found in 
EXPLORER-NAUTILUS vs the absolute value of the event 
time difference, obtained by subtracting the theoretical cu- 
mulative of accidentals from the cumulative of events found 
in proximity of the injection times. The common excitations 
are software injections of damped sine waves with / = 914Hz, 
t = lms. About 95% of the events with 4.2 < SNR < 4.6 
are retrieved with a coincidence window of ±25ms. 



FIG. 6: Responses of detectors to a damped sinusoid (h rss = 
1 • 10~ 18 Hz~ 1/2 , t = 30ms). The plot shows the amplitude 
response of Explorer and of Nautilus vs the wave center fre- 
quency. Markers indicate values measured via software in- 
jections, while the lines show the behavior predicted by an 
analytical model. 
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to damped sinusoids, linearly polarized and optimal ori- 
ented, with amplitude h rss = 1 • lO^^Hz^ 1 / 2 and decay 
time r = 30 ms, as a function of their central frequency 
/o- 

Analogously, figjT] shows the ratio of the responses 
(EXPLORER/NAUTILUS) for 4 different values of r: 
t is evident that, while this ratio is close to unity for 
short bursts, approaching delta pulses, it can vary up to 
a factor 3 for longer signals, since these two detectors, 
so similar under many aspects, do have narrow, not fully 
overlapping band widths. 

These considerations clearly show that the same exci- 
tation can produce different responses in the detectors. 
Since we are searching small, near threshold signals, we 
chose not to implement criteria that select events by re- 
quiring equal or comparable amplitude response, because 
this would only select short bursts, a priori rejecting 
longer signals for which the detectors are just as sen- 
sitive. Therefore, it does not appear sensible to apply 
a threshold based on the absolute H value of the candi- 
date events. This leads us to discard criterion C, that 
anyhow was adopted for calculating upper limits, not a 
goal for this paper. Both A and B criteria improve the 
efficiency to longer bursts with colored spectrum. Crite- 
rion B, instead, shows a better efficiency than A for short 
bursts, as this would give a higher SNR in the detector 
with better spectral sensitivity. Therefore in this search 
we adopted only criterion B. 
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FIG. 7: Responses of detectors to damped sinusoids with var- 
ious duration (r =1,3,10,30 ms). The plot shows the ratio of 
the responses (Amplitude Explorer/ Amplitude Nautilus) vs 
the wave center frequency. Markers and lines as in fig[6] 



C. Background estimation and fine tuning of the 
amplitude thresholds 

The FARs were estimated with the usual time-shifting 
technique, i.e. creating a large number of replicas of the 
same event lists with the event times of each list changed 
by offsets. In the past, this was done by choosing one 
shift step and producing replicas with offsets that are 
multiples of such step. 

The expected statistical distribution of accidental 
background can deviate from the Poisson distribution 
if clustering or non-stationarity affect the event lists. 
These deviations would go undetected if an uneffective 
time shifts is step chosen: large time shift steps do not 
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FIG. 8: Contour lines of the accidental rate (events/century) 
estimates in the AU-EX-NA (data stretch B)case. 



SNR threshold. 

The thresholds were tuned requiring the same SNR 
value for ALLEGRO, EXPLORER and NAUTILUS, and 
an higher SNR value for AURIGA, which has a noise level 
about 2 times smaller than the other detectors, thus bal- 
ancing the average sensitivities of the different detectors. 





Overlap 




SNR threshold 


Configuration 


time 




Stretch A 


Stretch B 




(days) 




2005-2006 




2007 




AL AU EX 


361.8 




4.3 


6.88 


4.3 


4.3 


6.24 


4.3 


AL AU NA 


390.6 




4.3 


7.1 


4.3 


4.3 


6.34 


4.3 


AL EX NA 


308.7 




4.94 


4.94 


4.94 


4.84 


4.84 


4.84 


AU EX NA 


301.9 




6.88 


4.2 


4.2 


6.14 


4.2 


4.2 



TABLE III: Time of simultaneous operation of the four three- 
fold configurations and SNR threshold for the various detec- 
tors . The different sets of values in Stretch A and Stretch B 
were chosen in order to yield the same FAR of 0.2/century in 
each search. 



correctly display possible correlation effects in nearby 
events, producing a poissonian behavior even if the ac- 
tual background of accidental coincidences shows a more 
complex structure at small time shifts. This effect prac- 
tically mimics a decorrelation of the data. On the other 
hand, a small shift step does not allow to probe a suf- 
ficiently large time lag, producing replicas that are not 
completely independent. 

In order to avoid this effect, we used, in each back- 



ground evaluation (see table IV) a series of 13 separate 
sets of evenly spaced time offsets, differing by the val- 
ues of shift steps {T off set } = {0.12,0.36, 0.60, 0.84, 1.08, 
1.32, 1.56, 1.80, 2.04, 2.28, 2.52, 2.76, 3.00} s. The corre- 
sponding maximum relative time shifts between any pair 
of detectors result in the range ±240 to ±6000 s. As an 
additional check, we verified the results of this procedure 
with other two alternative methods. The first one con- 
sists in assigning a random time offset, picked from a uni- 
form distribution to the event list of each detector. With 
this method, the distribution of time shift differences be- 
tween any pair of detectors is not uniformly distributed 
but concentrated around zero, thus better probing the 
zone around the zero-lag point. The second method is 
the one described in the appendix of [3] . In this case the 
accidentals are analytically estimated from the number 
of events in each single detector found in random sub- 
divisions of the overlap intervals. The results of both 
of these alternative methods agree, within statistical un- 
certainties, with those of the first technique, i.e. of the 
multiple time steps. 



1. Threefold coincidences analysis 

The target FAR of 0.2 events/century for each of the 
four different configurations of threefold coincidence was 
obtained by selecting the candidate events on the base of 



The selection was done starting from the lowest pos- 
sible SNR threshold, i.e. using all the exchanged data, 
and rising the threshold values till the target FAR was 
reached. Figj8] shows as an example how the FAR varies 
with the threshold values in AURIGA and in the pair 
EXPLORER-NAUTILUS for which a common value was 
adopted. This iterative process brought us to implement 
the SNR thresholds reported, for each configuration, in 
table HE 

We performed 12,006,000 time shifts per each of the 
13 steps cited above. 

All 13 different statistical distributions of the four con- 
figurations of three detectors are well fitted by a poisso- 
nian. The 52 mean values of the poissonian (i.e. the false 
alarms) are shown in table IV in the last row we report 



the mean and standard deviation per each detectors con- 
figuration. The experimentally determined ratio between 
standard deviation and mean is about 2.5 • 10~ 2 , larger 
than expected from a Poisson counting statistics, which 
would give 6 — 7 • 10 -3 . 

For each threefold configuration, we built the over- 
all distribution of the 13 ones, relative to different shift 
steps. The predicted probabilities of 0,1,2 and 3 acciden- 
tal counts per each detectors configuration are reported 
in table El 

The values in the first rows of each configuration are 
relative to the experimental background and were ob- 
tained calculating mean and standard deviation of the 
correspondent bins of 13 distributions each representing 
12,006,000 trials. The second rows of each configuration 
show the probability values in the hypothesis of poisso- 
nian distributions having means as in the last row of table 



IV The experimental and poissonian values are compat- 



ible within the uncertainties. 

We notice that to calculate the experimental values we 
did not consider that the same time shift can be obtained 
using different steps, but we verified that the effect of this 
overlap among the 13 used distributions is negligible. 



2. Fourfold coincidences analysis 

For the fourfold configuration, all the exchanged events 
of ALLEGRO, EXPLORER and NAUTILUS were used 
(thresholds as in table [TTJ) , while AURIGA data were se- 
lected so as to reach the desired 0.2 false alarms per cen- 
tury. AURIGA threshold resulted SNR = 4.96 for stretch 
A and SNR = 4.57 for stretch B (2007 data). As for the 
threefold case, also the analysis of the fourfold coinci- 
dence background was performed using the 13 time shift 
steps {T ff set }. The total number of independent trials 
for each shift was 11,094, 161. The last column of table 
IV reports the measured false alarms for each {T ff set } 
step. 

The last row of table [V] contains the experimental 
probabilities of obtaining 0,1,2,3 coincidences by chance. 
These values can be compared with those in the second 
line in the row, relative to a Poisson distribution of mean 
0.199 ± 0.003. 

Similarly to the threefold coincidence analysis, the ex- 
perimental probability is fully compatible with a poisso- 
nian distribution within the calculated errors. 



IV. RESULTS OF THE SEARCH AND 
CONCLUSIONS 

When the confidential time offsets were finally ex- 
changed, the true time analysis was performed and no 
real (i.e. true time) coincidence was found in any of the 
five parallel searches that, based on our background esti- 
mate, yielded a total FAR of 1 event per century. We also 
looked for coincidences on all exchanged data, i.e. with 
the amplitude thresholds of the data exchange shown 
in table [Ui| We found no quadruple coincidence and 20 
triples, well within the expected occurrence of accidentals 



(~ 22, as reported in table VI I . Since in this case the FAR 



is ~ 14/year, this analysis cannot identify single GW can- 
didates with a reasonable significance. In fact, we a pri- 
ori agreed that the only valid IGEC2 search for GWs was 
the former at FAR=l/century. Moreover, the total num- 
ber of coincidences found is well within the expectations. 
Nonetheless, as we had a priori decided, we make public 
[213] these coincidences to allow possible further analysis 
with the use of fresh data from other detectors or astro- 
physical observations. To allow a blind analysis of these 
data the published times of these coincidences have been 
shifted within ± 10 s, as performed internally by IGEC2. 
The signals with lowest amplitudes that IGEC2 can de- 
tect have an h rss of the order of a few 10~ 20 H z~ x / 2 , 
which is larger by more than one order of magnitude 
with respect to the current LIGO- Virgo network sensi- 
tivity for signals in the same frequency band [5j. Indeed, 
a simple calculation shows that, due to the long obser- 
vation time, IGEC2 can set an upper limit on incoming 
g.w. flux that is twice better than that set by LIGO S5 
run [HI IS], but only for amplitudes h rss > 10~ 19 ifz -1 / 2 . 

In this paper we have reported on a search for coinci- 
dent events on four resonant gravitational wave detectors, 
covering a data taking period of 17 month that produced 
no candidate event within our choice of confidence level. 
Since ALLEGRO ceased operating in 2007, this is a final 
report for the four antennas array, but the survey for a 
rare, highly energetic GW event does continue with the 
three detectors still in operation. Next searches will most 
likely be based on triggers from other observatories (neu- 
trinos, X or gamma rays, etc.) due to the lower statistical 
robustness of a network of three antennas, all located in 
the same geographic area. 

We believe that the IGEC2 survey can complement the 
observation time of the LIGO- Virgo network in the next 
years, until the long baseline interferometers resume their 
long term observations after the planned upgrades. 
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False Alarm Rates 




Shift 




(events / century) 




(s) 


AL AU EX 


AL AU NA 


AL EX NA 


AU EX NA 




AU EX NA AL 


0.120 


0.204 


0.234 


0.215 


0.203 




0.200 


0.360 


0.200 


0.196 


0.205 


0.195 




0.200 


0.600 


0.195 


0.199 


0.206 


0.196 




0.194 


0.840 


0.192 


0.199 


0.211 


0.205 




0.195 


1.080 


0.200 


0.200 


0.202 


0.202 




0.198 


1.320 


0.204 


0.199 


0.200 


0.207 




0.205 


1.560 


0.199 


0.198 


0.197 


0.203 




0.199 


1.800 


0.191 


n 1 qs 


n i qr 


n 9nn 




0.199 


2.040 


0.199 


0.197 


0.206 


0.206 




0.196 


9 9RH 


n i Q9 
u. lyz 


0.197 


0.198 


0.201 




n i qq 


2.520 


0.189 


0.206 


0.194 


0.203 




0.196 


2.760 


0.197 


0.203 


0.194 


0.210 




0.201 


3.000 


0.200 


0.201 


0.191 


0.204 




0.200 


mean, st.dev. 


0.197, 0.005 


0.202, 0.010 


0.201, 0.007 


0.203, 0.004 




0.199, 0.003 



TABLE IV: False Alarm Rates (FAR) of the four coincidences analyses (three threefold and one fourfold) , evaluated with about 
12 millions trials for each of the 13 time shifts steps of the set {T //set}: In each search the thresholds were chosen according 
to the values of table [TTTI 



Configuration 


P(N=0) 


P(N=1) 


P(N=2) 


P(N=3) 


AL AU EX 


0.998049 ±4.6- 10~ b 
0.998049 ±4.7- 10~ 5 


(1.949 ±0.046) ■ 10" 3 
(1.949 ±0.047) • 10 -3 


(2.02 ±0.49) ■ 10" b 
(1.904 ±0.091) • 10~ 6 


< 8- 10 _s 
(1.241 ±0.088) ■ 10~ 9 


AL AU NA 


0.99784 ±1.0- 10" 4 
0.99784 ±1.0- 10" 4 


(2.15 ±0.10) • 10" 3 
(2.19 ±0.10) • 10" 3 


(2.23 ±0.35) • 10 _e 
(2.34 ±0.22) • 10" 6 


< 8- 10 _s 
(1.69 ±0.24) ■ 10" 9 


AL EX NA 


0.998299 ±5.7- 10~ b 
0.998299 ±5.7- 10~ 5 


(1.700 ±0.057) • 10" 3 
(1.700 ± 0.057) • 10~ 3 


(1.49 ±0.29) • 10 _e 
(1.448 ±0.097) • 10~ 6 


< 8- 10" s 
(8.24 ±0.83) ■ 10 -10 


AU EX NA 


0.998325 ±3.4- 10"° 
0.998325 ±3.4- 10~ 5 


(1.674 ±0.034) • 10" 3 
(1.674 ±0.034) • 10~ 3 


(1.51 ±0.26) • 10 _b 
(1.40 ± 0.057) ■ 10~ 6 


< 8- 1(T S 
(7.85 ± 0.48) ■ 10~ 10 


AL AU EX NA 


0.998402 ±2.4- 10~° 
0.998403 ±2.4- 10~ 5 


(1.598 ±0.024) • 10" 3 
(1.595 ±0.023) • 10~ 3 


(2.1 ±3.9) • lO" 7 
(1.275 ±0.038) • 10~ 6 


< 9- 10"* 
(6.79 ±0.30) ■ 10 -10 



TABLE V: The first line of each configuration contains the experimental occurrence probabilities of 0,1,2 or 3 accidental 
coincidences calculated averaging the results of about 12 milions independent trials for each of the 13 shifts steps of the set 
{Toff set}- The second rows of each configuration show the poissonian probability values in the hypothesis of means as in the 
last row of table HVl 



Configur at ion 


Expected 
accidentals 
(events) 


Coincidences 
found 
(events) 


AL AU EX 


4.29 ±0.01 


3 


AL AU NA 


5.15 ±0.01 


5 


AL EX NA 


10.23 ±0.01 


8 


AU EX NA 


2.34 ±0.01 


4 


AL AU EX NA 


(7.66 ±0.01) • 10~ 3 






TABLE VI: Number of expected false alarms and of coinci- 
dences found in the search with minimum thresholds. 



